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On the Power-Law Tail in the Mass Function of Protostellar 
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1 INTRODUCTION 



, 1.1 IMF background 



Despite decades of detailed modeling of gravitational collapse and 
fragmentation, it is still not clear how certain physical conditions 
determine the mass of an individual star; a larger question is the 
origin of the distribution of stellar masses at birth, the stellar initial 
mass function (IMF). Additionally, the observed mass spectrum of 
protostellar condensations in molecular clouds appears to follow a 
power-law decrease, similar to that of stars, in the mass range above 
~ 0.5 M© (Motte, Andre, & Neri 1998; Testi & Sargent 1998; 
Johnstone et al. 2000). Hence, the fundamental shape of the IMF 
may be determined at the molecular fragment scale, even though 
the transformation of a condensation to star(s) is not a simple pro- 
cess. 

Starting with the expression for the Jeans (1928) mass, and 
working up to more sophisticated scenarios, there is a general ten- 
dency for multiplicative dependence on physical parameters (see, 
e.g., Adams & Fatuzzo 1996), leading to the claim that the cen- 
tral limit theorem should ensure that the IMF has a lognormal 
form. The formation of star-forming condensations in a molecu- 
lar cloud is surely affected by multiple effects of a turbulent flow 
field, and statistical variations of ambient magnetic field, density, 
pressure, and temperature. Numerical models of turbulent molecu- 
lar clouds typically yield a lognormal distribution of dense clump 
masses (e.g., Padoan, Nordlund, & Jones 1997; Ostriker, Stone, & 
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ABSTRACT 

We explore the idea that the power-law tail in the mass function of protostellar condensations 
and stars arises from the accretion of ambient cloud material on to a condensation, coupled 
with a nonuniform (exponential) distribution of accretion lifetimes. This model allows for 
the generation of power-law distributions in all star-forming regions, even if condensations 
start with a lognormal mass distribution, as may be expected from the central limit theorem, 
and supported by some recent numerical simulations of turbulent molecular clouds. For a 
condensation mass m with growth rate dm/dt oc m, an analytic three-parameter probability 
density function is derived; it resembles a lognormal at low mass and has a pure power-law 
high-mass tail. An approximate power-law tail is also expected for other growth laws, and we 
calculate the distribution for the plausible case dm/dt oc m 2 / 3 . Furthermore, any single time 
snapshot of the masses of condensations that are still accreting (and are of varying ages) also 
yields a distribution with a power-law tail similar to that of the IMF. 
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Gammie 2001; Klessen 2001). Although much additional physics 
remains to be included in these models, it may be reasonable to 
believe that the stochastic feature of a lognormal distribution will 
continue to hold. Hence, one can ask a simple set of questions about 
star formation: Is the IMF lognormal? If not, and an approximate 
power-law tail is inferred, then how is it generated? 

While Miller & Scalo (1979) did obtain a lognormal fit to the 
IMF, subsequent detailed studies find a region of power-law de- 
crease in the intermediate and high mass (M) range that is sim- 
ilar to the original result of Salpeter (1955), i.e., dN/d\nM oc 
M~ , where dN is the differential number of stars in the log- 
arithmic mass interval din M (see e.g., reviews by Scalo 1998; 
Meyer et al. 2000; Elmegreen 2001). Scalo (1998) emphasizes a 
varying power-law index in the intermediate and high mass range, 
with significant region to region variations reflecting intrinsic vari- 
ability and/or observational uncertainties. 



1.2 Lognormal or power-law distribution? 

The generation processes of lognormal and power-law distributions 
are often intimately connected. This phenomenon has been studied 
extensively in social science, biology, and computer science, e.g., 
in the study of power-law distributions of income (the Pareto Law), 
city sizes (Zipf's Law), and Internet file sizes. Useful models and 
discussion can be found in Reed (2002, 2003), Reed & Hughes 
(2002), Mitzenmacher (2003), and references within. It turns out 
that many seemingly small effects in the generative process of a 
lognormal distribution can lead to a power-law tail instead. Such 
an explanation seems particularly relevant for the IMF, since the 
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Figure 1. Estimated stellar masses in the ONC and best-fitting analytic 
curves. The squares represent the masses estimated by H97, binned in inter- 
vals A \og(M / Mq) = 0.2. The dashed line is the best-fitting lognormal 
distribution (eq. 0, with ^ = -1.40, a = 0.52) for M < 1 Mq, the 
dash-dotted line is the best-fitting lognormal distribution (p, = —1.08, a = 
1.16) for all masses. The solid line is the least squares best-fitting straight 
line for the region M 1 Mq ; it has slope F = d log N/ d log M 
■ 



power-law tail is typically measured to span even less than two 
decades of mass, hardly qualifying as evidence for a truly scale- 
free process. 

In this Letter, we propose a minimum-hypothesis explanation: 
even if the central limit theorem is relevant in determining the start- 
ing masses of condensations, the temporal effect of their subse- 
quent interaction with the environment skews the mass distribution 
into a power-law form. This paper is based on the statistical model 
summarized by Reed & Hughes (2002), and builds upon some el- 
ements of the model of Myers (2000). The basic idea is that the 
growth of masses by accretion coupled with a distribution of times 
for the growth yields a power-law tail in the mass function. It is a 
statistical approach, not dependent on the detailed values of phys- 
ical variables in different environments. While the mass of any in- 
dividual star is dependent on the particular initial values of various 
physical parameters, we are concerned here only with the shape 
of the mass distribution function; specifically the intermediate and 
high mass tail. 



2 LOGNORMAL FITS AND POWER-LAW TAIL IN THE 
IMF: A CASE STUDY 

In order to establish the minimum requirements for an IMF model, 
we review the data from the Orion Nebula Cluster (ONC) obtained 
by Hillenbrand (1997, hereafter H97). The ONC is a unique labo- 
ratory for the study of the IMF due to its relatively nearby location 
and direction away from the Galactic plane, allowing a deep sam- 
pling of the young stellar population in the vicinity of the massive 
star C. The large range of surface density from centre to ob- 
served edge of the cluster, and the large number of observed stars 
allows for useful statistics on the role of the environment (specifi- 



cally the density of stars) in the generation of the IMF. Furthermore, 
since most stars have estimated ages less than a few Myr (H97), 
they are likely observed near their birth sites. The IMF of H97 also 
shows a decline at lower masses, which we utilize to assess the 
viability of lognormal fits. 

Figure 1 shows a histogram of 696 stellar masses obtained 
by H97; these are the subset of observed stars which are reported 
to have a 70% or higher chance of membership in the ONC. We 
do not try to account for uncertainties in the mass determinations 
here, and simply seek a best-fitting lognormal probability density 
function 
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for the normalized masses m = M/Mq. For comparison with 
a sample of iVtot objects, the normalized number in each bin, 
dN/Ntot, equals f(m)dm — mf(m)dlnm, where dlnm is the 
width of a logarithmic mass bin. We conduct a parameter search of 
values of the mean fi and variance o 2 of In m. Minimization of the 
X 2 deviation of the observed values from the analytic curves yields 
best-fitting parameters n — —1.08 and a = 1.16. The dash-dotted 
line in Figure 1 is a lognormal distribution with these parameters. 
Clearly, the ability to fit the tail of the distribution is inconsistent 
with fitting the peak at the low mass end. The reduced \ 2 f° r this 
fit is x? = X 2 l v ~ 6-64 for u = 12 degrees of freedom, yield- 
ing a probability 4.75 x 10~ 12 that the data is drawn from the 
lognormal distribution. Alternatively, we may seek the best-fitting 
lognormal to a subset of the data; fitting only stars in bins centred 
at M < 1 Mq yields best parameters [i — —1.40 and o = 0.52. 
Figure 1 also shows this fit (dashed line), which has x? = 1-86 
and v = 3, yielding a probability 0.134 that the data is drawn 
from the distribution. A conclusion is that the IMF can be approx- 
imately fit by a lognormal at low mass, but has a power-law inter- 
mediate and high mass tail. The histogram yields a least squares 
best-fitting straight line of slope T = dlog iV/dlog M = —1.13 
for 1 Mq. 

Due to the large variation in stellar surface density in the ONC, 
we also look for systematic location-dependent effects on the IMF 
by dividing the stars into three radial zones, each containing about 
232 stars, enough for a meaningful statistical sample. Table 1 lists 
the values of x? from comparing the histogram of masses M < 
1 Mq in each zone with the best-fitting lognormal for all zones: 
H = —1.40, o — 0.52. The implied probabilities range from 0.287 
to 0.616. For M ^ 1 Mq, the best-fitting slopes to the power-law 
tail are shallower than —1 in radial zones 1 and 2, although there 
are significant uncertainties in the fits. Interestingly, if we exclude 
the (merely five) stars in bins centred at M > 10 Mq, all slopes 
become somewhat steeper than —1. Four of the five excluded stars 
are in zone 1 (these are the Trapezium stars) and one is in zone 2. 

The ONC data implies a generation of essentially the same 
power-law in all regions, despite large differences in stellar density 
(see Table 1). The five most massive stars, four of which are very 
near the core of the ONC, provide the only significant evidence for 
a region-dependent IMF. Their presence and location may be due to 
events occurring after star formation, since the local relaxation time 
decreases toward the cluster centre due to increasing density. How- 
ever, the simulations of dynamical mass segregation by Bonnell & 
Davies (1998) imply that they could not have formed very far from 
the cluster centre. In any case, we bypass the nagging questions 
about the origin of these five stars and focus instead on an expla- 
nation for the generation of a universal power-law tail in the mass 
function of all other stars. 
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Table 1. Summary of best-fitting parameters. 



Zone 


Surface Density (pc 2 ) 


M < 


1M Q 


M ^ 1M Q 


1M < M ^ 10 M 


lognormal x 2 


probability 


best slope 


best slope 


all (0 < r < 17.70 pc) 


0.71 


1.86 


0.134 


-1.13±0.15 


-1.53 ±0.27 


1 (0 «C r < 4.32 pc) 


3.94 


0.598 


0.616 


-0.85 ±0.21 


-1.36 ± 0.39 


2 (4.32 < r < 8.65 pc) 


1.32 


1.25 


0.287 


-0.89 ±0.16 


-1.12 ± 0.22 


3 (8.65 < r < 17.70 pc) 


0.31 


1.17 


0.318 


-1.10 ±0.45 


-1.10 ±0.45 




Figure 2. Evolution of a distribution of normalized masses m undergoing 
accretion growth dm/dt oc m (eq. |2|). The dashed line is an example ini- 
tial lognormal distribution with fio = — 1.40, cro = 0.52, and the circles 
represent a histogram (binned with A log m = 0.2) of a random sample of 
10 6 objects drawn from the distribution. The dash-dotted line is the analytic 
expression for the distribution (eq. Q with /t = /to + 1) if the masses grow 
for a fixed time t = 7 . The solid line is the analytic density function 
(eq. Q with fj,o = —1.40, (To = 0.52, a = 1) if accretion lifetimes have 
density function /(t) = Se~ st , and 5 = 7. The diamonds and squares rep- 
resent the histogram of the masses in the original distribution if they grow 
for a fixed time or with an exponential distribution of times, respectively; 
they necessarily overlap the analytically predicted curves. 



3 A STATISTICAL MODEL OF GROWTH BY 
ACCRETION 

3.1 Basic model 

According to the central limit theorem of statistics, if the mass of a 
protostellar condensation M c = fi X f% X ... X /jv, then the distri- 
bution of Af c tends to a lognormal regardless of the distributions of 
the individual physical parameters fi (i = 1, N), if N is large. 
Depending on the specific distributions of the fi, a convergence to 
a lognormal may even occur for moderate N (Ioka & Nakamura 
2002). However, the stellar IMF results from further evolution be- 
yond the condensation phase, so that the mass of a star may be 
characterized as M s = M c x /f rag x / ou t, where /f rag is the frac- 
tion of the condensation mass that goes into a star after any possible 
fragmentation process, and / ou t is the fraction of mass that remains 
after mass loss due to a protostellar outflow. The multiplicative na- 



ture of these processes seems consistent with theoretical models 
(Shu et al. 1999; Meyer et al. 2000), and imply that the distribution 
of M s should be lognormal if that of M c is lognormal as well. 

Here we explore the idea that the power-law tail is due to the 
interaction of a condensation with its environment in the time be- 
fore it collapses down to stellar dimensions. The environment of a 
condensation is certainly a large reservoir of mass, and there may 
be dynamical influences due to the passage of ambient pressure 
fluctuations, including shocks. If the accretion is a continuous time 
multiplicative process, 

dm 
~dt 

where 7 is a growth rate and mo is an initial mass. If the values 
of mo are drawn from a lognormal distribution with mean fio and 
variance <r 2 , the distribution of masses f(m) at any later time t will 
still be described by equation 0, but with 



= 7m =>■ m(t) = mo exp(7 1), 



(2) 
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However, all condensations may not accrete for the same time 
period, and a distribution of accretion times t will in fact skew the 
final distribution away from a lognormal. A simple example of a 
distribution of times is if the probability of stopping accretion is 
constant in time. A constant "death" rate S for accretion leads to an 
exponential distribution, with probability density function f(t) = 
Se~ st . In this case, the probability density function for masses is 
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is the complementary error function, we find 
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where a — 5/7 is the dimensionless ratio of "death" rate to 
"growth" rate of condensations. Equation represents a new 
three-parameter (/xo , 00 , a) probability density function which 
tends to the form m _1_a for large masses, but is modulated by 
the complementary error function (which goes to for very low m 
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and to 2 for large m) so that it tends toward zero at very low masses. 
It is therefore similar to a lognormal distribution with a power-law 
tail. (Note that it reduces to a pure power-law distribution if all stars 
start with the same mass mi, i.e., mo is described by a delta func- 
tion, so that fj,o = In mi and ao — in eq. Q) We believe this 
function can be used to fruitfully model the IMF. 

Physically, we expect a ~ 1 since both 5 and 7 are rates 
controlled by the external medium and are expected to be approx- 
imately the inverse external dynamical time t^ 1 = (Gp cx t) , 
in which p cxt is the mean density outside the condensation. My- 
ers (2000) has derived an example of exponential growth with 
7 w tT 1 , assuming geometric accretion (see eq. (5)) and that every 
condensation increases its internal turbulence continuously (as the 
mass increases) exactly as needed to remain in a critical Bonnor- 
Ebert state. 

Figure 2 shows a sample initial lognormal distribution (jj,o = 
— 1.40, (To = 0.52), and a distribution at a later time under the 
assumptions (1) multiplicative growth (eq. (2)) in which all con- 
densations accrete for a fixed time t = 7 -1 , and (2) multiplicative 
growth with an exponential distribution of times f(t) = 8e~ st . 
Under assumption (1), the new distribution is again a lognormal 
with fA = fxo + 1> but in case (2) it is given by equation and has 
a power-law tail with index a = 1 if 5 = 7. 

Recently, Reed (2002, 2003) has derived a probability den- 
sity function that is even more general than equation 0. It as- 
sumes the growth law known as Geometric Brownian Motion: 
dm/m = 7 dt +<7i dw, where dw is a "white-noise" random num- 
ber drawn uniformly from the range [0, 1]. In this case, an initially 
lognormal distribution characterized by /10 and ao also remains 
lognormal after a fixed time i, but with variance a 2 — (Tq + <j\ t, 
and mean fj, = fj,o + (7 — cr?/2) t. For an exponential distribution 
of lifetimes, Reed derives an analytic four-parameter density func- 
tion, more complicated than equation Q, but which also has the 
asymptotic dependence /(to) oc m -1 ~ a for large m. We refer the 
reader to the above-mentioned papers for details. 

While an exponential distribution of accretion termination 
times guarantees a power-law tail in the mass distribution, it is 
also true that all condensations need not start accreting at the same 
time. If the "birth" of condensations is distributed uniformly in time 
throughout the life of a cloud, and termination times obey an expo- 
nential distribution, then the number of "living" (i.e., still accreting) 
condensations at a fixed time t of observation will also obey an ex- 
ponential distribution of times r since birth, i.e., f(r) = Se~ Sr . 
Furthermore, if the "birth" of condensations is also increasing ex- 
ponentially in time with a growth rate j3, then the distribution of 
times r since birth will be /(r) = (/3 + 5)e~" 3+ *' T . Hence, ob- 
servations of accreting condensations made at a single snapshot 
in time, e.g., Motte et al. (1998), should yield a power-law index 
a! = q + /3/7, so that the slope is equal to that of the IMF if = 
(or j3 <C 7), but steeper otherwise. 



3.2 A variation on the basic model 

Accretional growth proportional to the instantaneous condensation 
mass is a plausible but unproven assumption. Myers (2000) used 
such a relation for the growth of condensations, but had to assume 
a very specific rate of growth of internal turbulence as accretion 
proceeded. 

Here, we look at an alternative scenario, also based on 
Bonnor-Ebert spheres, but assuming that the internal velocity dis- 
persion remains fixed during accretion. It is instructive to plot the 
radius versus mass (R — M) relation for the sequence of equilib- 
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Figure 3. Radius (R/Ro) versus mass (M/Mq) relations for a sequence 
of Bonnor-Ebert equilibrium states for constant external pressure P e xt 
(dashed line). Here, R = ^/(G^P^ 2 ) and M = cf/(G 3 ^ 2 P^). 
The plotted sequence extends somewhat beyond the critical mass M CT - lt = 
1.18 Mo. The dash-dotted line is the corresponding sequence for an exter- 
nal pressure 3 P C xt ■ The solid line with arrows represents a possible evolu- 
tionary sequence of a condensation during mass accretion, terminating with 
collapse as the external pressure suddenly rises to a value incompatible with 
an equilibrium state. 
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Figure 4. Evolution of a distribution of normalized masses m undergoing 
accretion growth dm/dt oc m 2 ' 3 (eq. HOH . The dashed line is an exam- 
ple initial lognormal distribution with no = — 1.40,<ro = 0.52. A sam- 
ple of 10 6 masses drawn from this distribution are binned in increments 
A log m = 0.2 and represented by circles. The diamonds (connected by a 
dash-dotted line) are a histogram of masses after they have all accreted for 
a time t = 7^ 1 , and the squares (connected by a solid line) are a histogram 
of masses if the accretion lifetimes have density function f(t) = 5e~ st , in 
which 8 = 71 . 
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rium states of increasing mass, but with fixed internal sound speed 
c s and external pressure P oxt . For small clouds, the R — M relation 
is found analytically from the Bonnor-Ebert theory (Bonnor 1956; 
Ebert 1957) to be 



R 



3ci 



47rP cxt 



1/3 



M 



1/3 



(8) 



However, this approximation also yields radii which are within 
20% of the actual value for all M/M < 1, where M = 
4/(G 3/2 P e L x / t 2 ) is the natural unit of mass in the Bonnor-Ebert 
problem. Figure 3 shows the Bonnor-Ebert R — M sequence as 
well as the corresponding relationship for a higher external pres- 
sure 3P cxt . 

In this simplified model, condensations start out somewhere 
on the equilibrium R — M sequence, increase their size R accord- 
ing to the curve as M grows by accretion, then make a jump down- 
ward after an ambient pressure fluctuation, as shown schematically 
by the solid line in Figure 3. After the downward jump, no accessi- 
ble equilibrium state may be present, and the condensation begins 
dynamical collapse. 

If R exceeds the Bondi (1952) radius Rb = GM/c SjCXt , 
where c s , cx t is the sound speed in the external medium, then a sta- 
tionary condensation undergoes geometric accretion 



dM A R 2 

— — = 4ttR p cxt Cs,ext- 

Equations JSJ and {5J then lead to 

2/3 -> m(t) 



dm 
~~dt 



71 m 



™V3 i Ti + 



(9) 
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w here m = .1 / : .V/h is a normalized mass, mo is the initial normal- 



ized mass, and 71 = (367T) 1 ' 3 t d 1 = 4.84 t d 1 is a growth rate. 



Figure 4 shows the evolution of the probability density func- 
tion f(m) from an initial lognormal distribution (same as in Fig. 
2) to a later distribution under the assumptions (1) all conden- 
sations grow for a fixed time t = 7-j" , and (2) condensations 
grow for times that are distributed exponentially, with "death" rate 
8 = 71. The first case leads to a profile that resembles a lognor- 
mal, but is narrower than the initial distribution, and the second 
case leads to a distinctly non-lognormal profile which may be de- 
scribed as a broken power-law with two distinct regions; best-fitting 
slopes r = d log N/d log m are equal to —1.0 in the interme- 
diate mass range (1 ^ m < 10), —2.4 in the high mass range 
(10 ^ m < 100), and —1.6 for all m ^ 1. Interestingly, a bro- 
ken power-law with a steeper slope at high masses is reminiscent of 
some descriptions of the field-star IMF (Kroupa & Weidner 2003). 



4 CONCLUSIONS 

We believe that the basic features of this model will be preserved in 
more complicated and realistic scenarios; that is, even a lognormal 
distribution of condensation masses at birth coupled with subse- 
quent accretion growth (proportional to some power of mass) in 
which the accretion time obeys its own distribution (approximately 
exponential), will result in a power-law tail. 

Zinnecker (1982) has presented another model of accretion 
growth in which a distribution of initial masses mo undergo Bondi 
accretion (m oc m 2 ); the exponent greater than 1 means that small 
differences in the mo values are magnified over time. After a fixed 
accretion time for all objects, when m 2> mo, /(m) tends to a 
power-law form with index —2. For comparison, the mechanism 



of exponential distribution of lifetimes studied in this paper yields 
a power-law tail even if the initial masses are all the same. The 
general combination of accretion growth (due to the Bondi pro- 
cess or geometric accretion) and a distribution of accretion life- 
times (where both processes are controlled by the dynamical time 
of the parent system) may play a role in determining the IMF in 
detailed global simulations of star formation (e.g., Bate, Bonnell, 
& Bromm 2003) in which stars are formed throughout the dura- 
tion of the simulation. Other ideas related to the spatial structure 
of turbulence may also contribute to power-law fragment masses, 
e.g., thresholded regions in density fields with lognormal probabil- 
ity density functions may have power-law clump mass distributions 
(Pudritz 2002; Elmegreen 2002). However, we believe that the tem- 
poral effect studied in this paper provides a compelling part of the 
explanation for the power-law tail in the IMF. 
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